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Abstract 



I explain a generalization of Bjorken flow where the medium has finite transverse size 
and expands both radially and along the beam axis. If one assumes that the equations of 
viscous hydrodynamics can be used, with p = e/3 and zero bulk viscosity, then the flow I 
describe can be developed into an exact solution of the relativistic Navier-Stokes equations. 
The local four-velocity in the flow is entirely determined by the assumption of symmetry 
under a subgroup of the conformal group. 
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1 Introduction 

Bjorken flow [1J is an attempt to describe the average motion of partons after a collision 
of heavy ions. Three assumptions that go into the treatment of [TJ are approximate boost- 
invariance along the beamline near mid-rapidity, translation invariance in the transverse 
plane, and rotation invariance in that plane. If Minkowski space is parametrized^ as 

ds 2 = -dr 2 + r 2 drf + dx\ + x 2 ± d<p 2 , (1) 

then translation and rotation invariance in the transverse plane imply that nothing can 
depend on x± or <p, while boost-invariance (if exact) implies that nothing depends on rj. 
Together with invariance under reflections r\ — > —rj, these symmetries imply that the local 
four-velocity vector is = (1,0,0,0) in the (r, r],x±,4>) coordinate system. One does not 
need any information about the equation of state to determine the four-velocity profile. 
Symmetry considerations alone fix it. The medium need not even be equilibrated for w M = 

1 More properly, the coordinates |lj cover a wedge of R 3:1 which is the causal future of the collision plane. 
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(1,0,0,0) to be a well-motivated choice for the average four-velocity in the mid-rapidity 
region. If the medium is equilibrated and has equation of state p = e/3, then the energy 
density in the local rest frame scales as l/r 4//3 . 

Translation invariance in the transverse plane is, of course, highly unrealistic since the 
nucleus is only about 13 fm across. Treatments based on the Bjorken picture often assume 
that the medium has no average local velocity in the x± direction (radial flow) until after 
it thermalizes. This is probably wrong, even for perfectly central collisions, and one might 
worry that it significantly distorts the subsequent hydrodynamical flows on which much of 
heavy-ion phenomenology depends. Indications that the absence of initial radial flow could 
be problematic can be found, for example, in (2J [3] ; see also [U El E] • The question naturally 
arises, can we estimate in some way deviations from the Bjorken picture, based on the finite 
size of the colliding nuclei, which lead to non-zero u x± 7 The aim of this paper is to present 
a generalization of Bjorken flow which does exactly that. 

The approach I will follow is based mostly on symmetry considerations, and it requires 
that the initial state and the dynamics of the collision perfectly respect relativistic conformal 
invariance. Also, I am limited to perfectly central collisions. Obviously, these assumptions 
are by no means exactly satisfied in real collisions. However, Quantum Chromodynamics 
(QCD) processes well above the confinement scale are close to conformally invariant because 
the coupling runs only logarithmically with scale. So conformally invariant dynamics, es- 
pecially in early stages of the collisions, is an interesting start point. The key input which 
allows me to extract a velocity profile solely from symmetry considerations is the assump- 
tion that translations in the transverse plane can be replaced by certain special conformal 
transformations which, in combination with rotations around the beamline, form an 5*0(3) 
subgroup of the full conformal group. The significance of this £0(3) subgroup was previ- 
ously noted in calculations [7J based on the gauge-string duality [HI El EE]- Indeed, ideas in 
[71 ITT] (see also [121 H3]) strongly underlie the proposal in this paper. But I do not depend on 
any of the dynamical information that the gauge-string duality provides. This is in contrast 
to [14] , where a gauge-string dual description of Bjorken flow was found by constructing an 
approximate solution to Einstein's equations in AdS^. The main calculations in this paper 
rely only on symmetries and hydrodynamics. 

The organization of the rest of the paper is as follows. In section [2] I will explain the 
aspects of conformal symmetry that I am going to use. In sectional use conformal symmetry 
to pick out a special four-velocity profile. In section [4] I will use conformal symmetry again 
to help find an energy density which, along with the special four-velocity profile, solves 
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the equations of inviscid, conformally invariant relativistic hydrodynamics. In section [5] I 
extend the result to the case of non-zero shear viscosity. The exact solution of the Navier- 
Stokes equations that I find generalizes the result e oc 1/r 4 / 3 for Bjorken flow. Unless the 
viscosity is identically zero, there is a pathology in the solution at very early times and/or 
very large radius. This pathology can be understood as the breakdown of a hydrodynamical 
description. In section [6] I characterize the colliding shock waves which respect the same 
conformal symmetries as the special four-velocity profile after the collision. In section [7] I 
use results on total multiplicity to estimate numerical parameters in the hydrodynamical 
flows found in previous subsections. In section [8] I explain how some of the features of 
the flow are more transparent when one maps it to 5 3 x R, which is the covering space of 
Minkowski space, R 3 ' 1 . Readers wishing to skip the motivations and technical detail can 
find a brief summary of the main results in section |9j 



2 An 50(3) subgroup of the conformal group 

The conformal group in four dimensions is an extension of the Poincare group 750(3, 1) 
to 5*0(4,2). The generators of the Lie algebra underlying 5*0(4,2) are (essentially in the 
notation of [15]): 

• Translations: £ M = d 1 where d 1 is constant. 

• Rotations and boosts: £ M = u^ u x u where is constant and anti-symmetric. 

• Scale transformations: £ M = x^. 

• Special conformal transformations: £ M = x u x u b fl — 2b u x u x fl where is constant. 

Obviously, the first two classes of symmetries belong to the Poincare group, 750(3, 1). 
Geometrically, they are Killing vectors, satisfying C^g^ u = where C$ denotes the Lie 
derivative and is the standard metric of R 3 ' 1 with mostly plus signature^] The vector 
fields associated with scale transformations and special conformal transformations are not 
Killing vectors, but rather conformal Killing vectors, satisfying 

^g^ = \(y x e)g»v (3) 

2 Readers unfamiliar with Lie derivatives can understand the main points by interpreting C% as an operator 
which measures how a tensor changes due to a coordinate transformation — > af^+^j treated to linear order 
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Often, I will follow conventions of differential geometry in specifying a vector in terms of a 
combination of derivatives: for example, £ = a^gfjr corresponds to £ M = a M . 
As reviewed in the introduction, the symmetries of Bjorken flow are: 

1. Rotation invariance around the beamline, £ = x 1 -^ — ^gfr- Of course, this is an 
exact symmetry only in the vanishingly rare case of a perfectly head-on collision. 

2. Translation invariance in the transverse plane, £ = and £ = J^. This is the 
symmetry I am most interested in relaxing, since it forbids radial flow in the transverse 
plane. 

3. Boost invariance along the beamline, £ = x 3 t| + t-£%. This symmetry is the key feature 
of Bjorken's treatment. It is based on a model where the Lorentz-flattened nuclei largely 
pass through one another and leave behind a medium in the wedge |x 3 | < t. The boost 
invariance is supposed to hold only in some neighborhood of mid-rapidity — which is to 
say, for \x 3 /t\ not too close to I. 

In order to show that these three symmetries, together with symmetry under the Z 2 group 
reflecting x 3 — > —x 3 , completely fix the local four- velocity u^, the equations we solve are 
[£, u] = where £ is any of the symmetries in the numbered list above. The only solution 
consistent with the Z 2 symmetry is u = J^, where we pass to coordinates (r, 77, x±, 0) defined 

as 

r = \ft 2 — (x 3 ) 2 i] = arctanh 



x 3 



x± = a/ (rr 1 ) 2 + (x 2 ) 2 (f) = arctan — 



^ < 4 > 

x 1 



A formal phrasing of what makes u = special is that it is the only timelike unit vector 
which is invariant under the I SO (2) symmetry of the transverse x^-x 1 plane, the 5*0(1, 1) 

in The main properties of are linearity, the Leibniz rule for products, and the following definitions: 

C^4> = ^ -^—^ for scalars <p 

Ctv? = £ A 7^r - uSrr- for vectors t/ 4 (2) 

ax A ox A 

% ^ A & +CA fS forl-formsO, 

While is usually defined in terms of partial derivatives d/dx^, one may instead use the covariant derivative 
V/t, defined in terms of the Christoffel connection, without changing the results. A particularly useful result 
is C^g^u — V^i, + Vj,£ M for the metric tensor g^ u . 

It is also good to note that = —C u ^. One often denotes C^u by [£, u], also known as the Lie bracket, 

or simply as the commutator of £ and u. 
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boost symmetry in the x 3 direction, and the Z 2 symmetry that reflects through the collision 
planej^] 

The ISO(2) x 5*0(1, 1) symmetry also implies that e, the energy density in the local rest 
frame of the medium, can depend only on r. But in order to actually determine the functional 
dependence e(r), one needs much more information. For example, if hydrodynamics is 
assumed to be valid, then one needs conservation of the stress-energy tensor, constitutive 
relations, the equation of state, and a choice of shear and bulk viscosity. 

I am interested in modifying the symmetry constraints in order to accommodate finite 
transverse size and radial flow. After some thought starting from the treatment [7J [TTJ of 
colliding shocks in AdS'sJ^] I concluded that the way to proceed is to replace translation 
invariance by invariance under 



9 2 
OX 



'id d 

Ox^ ox 1 



(5) 



where i = 1 or 2 and q is a quantity with dimensions of inverse length. The vector £ 14 (and 
analogously £24) is so named because it is inherited from a symmetry of AdS$, embedded in 
R 4 ' 2 by the equatior^] 

-(X- 1 ) 2 - (X ) 2 + (X 1 ) 2 + (X 2 ) 2 + (X 3 ) 2 + (X 4 ) 2 = -1 , (6) 

and this symmetry acts on X 1 and X 4 as an ordinary rotation. Even without reference to 
AdSs, ([5j is a more obvious thing to try than any other conformal transformation outside 
of IS0(3, 1), because it is a combination, for % — 1, of translation by a M = 51 with a special 
conformal transformation with 6 M oc It is easily checked that the commutator of £14 
and £24 is proportional to Jr, which is rotation around the beam axis. Indeed, £14, £24, and 
generate an 50(3) subgroup of 50(4,2) which commutes with the 50(1,1) subgroup 
generated by £ = x 3 -^ + tjj^§. The particular choice of the 50(3) is controlled by the 
parameter q, and as q — > 0, the 50(3) degenerates to the 750(2) symmetry of Bjorken flow. 
Let's call the group generated by ^4 and £^ 50(3) g to remind ourselves that it depends on 



3 Relaxing our insistence on the Z 2 symmetry while preserving the continuous symmetries would allow 
non-zero vP . This is a mild generalization because one could perform an overall boost to restore the Z2 
symmetry. 

4 Already in 111) the hope was expressed that heavy-ion collisions might be approximately invariant 
under an 0(3) symmetry (or an 0(2) subgroup for off-center collisions) generated as explained here. Failure 
of motivation and insight, in some combination, prevented me from actually writing down the key equation 



(131 for almost two years. 

J More properly, AdS^ is the universal cover of the solution space to (|6 
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q in an interesting way. 

In order to consider S0(3) q symmetry further, it will be convenient to recast one of the 
conformal Killing vectors ^ in terms of the coordinates Q: 

r - c 02 ,< 9 ,22,22n d 1 + q 2 T 2 - q 2 x\ d . . 

C = 4i4 = 2g rx± cos 0— + (1 + q r + q x ± ) cos 0- sin 0— . (7) 

3 A special four-velocity profile 

A definite problem is now almost formulated: I want to find a four-velocity profile in R 3 ' 1 
which respects SO(3) q x 5*0(1, 1) x Z 2 symmetry instead of the ISO(2) x S0(1, 1) x Z 2 
symmetry of Bjorken flow. (As before, the Z 2 symmetry is the reflection x 3 — > —x 3 , or 
equivalently rj — y —rj.) I expect that symmetry constraints will completely determine w M . 
Invariance under J^, J^, and rj — > —rj requires 



u = cosh k(t, x±)-^- + sinh k(t, x_|_) 



<9t <9x_i_ 

Readers familiar with the literature on hydrodynamical treatments of heavy-ion collisions 
will recognize ^ as a rewriting of the standard parametrization 

f d d d\ 



dr x dx v dy 



where ^ 

7± = y^=^= and u ± = ^v 2 x + v 2 . (10) 

For clarity of comparison with the literature, in ^ and (10) I have parametrized the trans- 
verse plane with (x,y) instead of (x 1 ,^ 2 ). The quantity v± is known as the transverse 
velocity. Comparing ^ and ^ leads immediately to 

u ± 

u T = cosIik = 7j_ — = tanhfi; = vj_ . (11) 

// 

Clearly, the vector u in (8) commutes with and J^: this is what "invariance under 
and means. But there is no choice of k(t, x±) that will get u to commute with (. The 
closest one can come is to choose 

2q 2 rx± 

K = arctanh — — j-^— — , (12) 
1 + q 2 r 2 + q z x\ 
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in which case one finds 



£ cMm = ^(V a C a K. (13) 



The result (13) at first seems discouraging, but in fact it is ideal. Let's refer to a tensor 
Q V uxu2 -- as having ^-weight equal to a if 

W^ = -t(VaC a )Q^- (14) 



According to (13), w M has weight —1, and according to ([3]), the metric has weight 
—2. Thus the projection tensor g fJiU + u fJi u l , has ^-weight —2. This projection tensor enters so 
ubiquitously into hydrodynamical equations that we should be delighted to see it transform 
as simply as possible under the symmetry generated by (. As I will explain in section [8j 
£- weight is closely related to the more general notion of conformal weight. 



4 The inviscid case 

Having decided upon a four-velocity profile, the next step is to consider what the energy 
density e should be. Let's start by requiring that e is invariant under and J^: thus 
e = e(r, x±). In the case q — 0, demanding that also (^d^e = would lead immediately 
to the conclusion that e depends only on r. It is tempting to require C^dne = even when 



q 7^ 0, but the experience of finding (13) instead of the more obvious condition C^u^ = 
is a hint that a slightly more elastic notion of symmetry is appropriate. Let us therefore 
consider the equation 



Ae = -f(V A C A )e- (15) 
Understanding e as a scalar quantity, this equation can be read as saying that the function 



e(r, x±) has (^-weight a. It is easy to see that the general solution to (15) is 



e= i M where g = 1 ~ gV + ^ (16) 
T a y 2qr V ' 

and e(g) is an arbitrary function. If one works to leading order in small q, one has g oc 1/r, 
and ( Jl6| reduces to the original conclusion that the energy density should depend only on 

T. 

Some dynamical information is required to pin down what e(g) is. The standard equations 
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of viscous relativistic hydrodynamics are 



(17) 



where 



_ p ap 



v a c 



To complete the equations, we must specify how the pressure p, the shear viscosity 77, and 
the bulk viscosity £ depend on e. Let's start with the simplest case: 



V 



v = C = o 



(19) 



The stress tensor is then traceless, as conformal invariance demands. The conservation 



equations (17) overconstrain e because both the v = t and v = x± equations are non-trivial. 



By inspection I found that these equations are consistent iff a = 4jj Then they reduce to 

dloge 8(7/3 



dg 



whose solution is 



(20) 



(21) 



Using (16), one finds immediately 



(2g) 



8/3 



^ 4/3 [1 + 2g 2 (r 2 + x\) + g 4 (r 2 



x 



2 \2 



4/3 ' 



(22) 



where eo is an integration constant. Recall that through mn and (12), we have already 



completely fixed u^. Also note that if we take q — > with eo<7 8//3 held fixed, then from (22) 



we recover the standard result e oc I/Y 4 / 3 for Bjorken flow. Thus (22) together with (18 



and (12) provide an exact solution of the equations of relativistic inviscid hydrodynamics 
with p = e/3 which becomes Bjorken flow in the limit q — > 0, but at finite q describes a 
medium which has integrable falloff in the x± direction. It is amusing to see an additional Z 2 
symmetry emerge which seems to have nothing to do with conformal symmetry: Based on 



6 In retrospect it seems obvious that a — 4: then if e is dimensionless, ( 22 ) leads correctly to the conclusion 
that e has dimension 4. 
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8D, (12), and (22), the quantities r 4//3 e, k, u t , and u 1 - are all invariant under the exchange 



of r and x±. 



5 The viscous case 

The obvious thing to try next is to generalize to non-zero viscosity while preserving conformal 
invariance of the theory: that is, 



V 



V 



H e 3 / 4 C = 0. 



(23) 



The e 3 / 4 dependence is necessary because only then is H dimensionless. In searching for a 
solution with non-zero shear viscosity, I'm not going to change w M at all: recall that is fixed 
entirely by the requirement that u M should have ^-weight equal to —1. I'm also not going to 
change the requirement that e should have ^-weight equal to 4: this is just the generalization 
to non-zero q of the condition that e should be constant across the transverse plane. In other 



words, I'm going to assume my, (12), and (22) with a = 4, and plug everything into (23) 



and (17) to get equations for e(g). Happily, the conservation equations are still consistent 
with one another for non-zero Ho- The equation they imply for i(g) is simpler to state in 
terms of 

T(g) = 0{g), (24) 
which is related to the local temperature of the fluid. The conservation equations imply 



,dT 



3(1 + g 2 ? 12 ^ + 2gV^+?T + g 2 R = 



(25) 



whose general solution is 



T(g) 



To 



+ 



1-(1+/) 1/6 2^1 



1 1 3 

2' 6' 2' 



(26) 



where To is an integration constant and 2-F1 denotes a hypergeometric function: 

2 F 1 (a,/3;r,z) 



aft a(a + l)/3((3 + 1) z 2 
c 7(7 + 1) 2 



q(q + l)(a + 2)/3(/3 + l)(/3 + 2) z 3 

7(7+1) (7 + 2) 3! + " 



(27) 
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There is a pathology in the solution (26): T is negative for large enough g. In fact, t has a 



single real root g*, and T < when g > g*. In other words, the T < pathology arises when 
t is too small and/or x± is too big. It shouldn't dismay us unduly since a similar pathology 



already arises in the q — > limit. To see this, expand (26) at large g (which is equivalent to 
small q) to find 



T(g) 



T — H - 



r( i/2)r(-i/3) 

2T(l/6) 



(2g) 2 / 3 H c 



T 



i/3 2r 



+ 0(g~ 



(28) 



Scaling q — > and Tq — >• oo so that Tog 2//3 remains finite, one finds 



e t 



T 

T 



rV3 



Ho 
2r 



(29) 



It can be checked that all the corrections to (29) from the terms labeled 0(g ) in (28) come 



with positive powers of q after the scaling just mentioned. Therefore, (29) combined with 
the q — y limit of the four-velocity profile, namely u = d T , forms an exact solution to the 
equations of conformal viscous hydrodynamics. This is viscous Bjorken flow. The pathology 
at r = t*, where 

/U_ \ 3/2 

(30) 



/Ho 
V2e 



can be understood as an indication that as r approaches r* from above, eventually hydro- 
dynamics cannot be used, because the shear is so strong that the viscous correction is more 
important than the pressure. In such a situation, it must be expected that higher derivative 
corrections also become important. The difficulties one anticipates as g approaches g»* from 
below in the finite q case are essentially the same as when r approaches r* from above for 
q = 0. 



6 Before the collision 

I'd now like to inquire what kind of non-hydrodynamical initial state might lead, at least in 
some approximation, to the hydrodynamical flow that I explained in the previous subsections. 
A natural ansatz in heavy-ion collisions is to assume that the state before the collision can 
be described as sum of left- and right-moving parts, each of which moves at the speed of 
light: 

T U u T uu (u : 3?j_) T vv T vv (v , x± ) , (31) 
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where 

u = t — x 



' v = t + x 3 . (32) 



All components of the stress tensor other than T uu and T vv are assumed to be zero in the 
coordinate system (u,v,x±,<f>) on R 3,1 . 

In order for the initial state (31) to lead approximately to an 50(3) g x SO(l, 1) x Z 2 



invariant hydrodynamical flow, it would help for the initial state to preserve as much of 
this symmetry group as possible. The 50(1, 1) boost symmetry can't possibly be preserved, 
because boosting a lightlike collision changes the total energy of the two participants mul- 
tiplicatively while preserving the product of their energies. 50(1, 1) symmetry is supposed 
to emerge in the mid-rapidity region through some post-collision dynamics. But 50(3) g 
symmetry, as the analog of 750(2) symmetry in the transverse plane, is something we 
might sensibly demand of the initial state in a perfectly central collision. Let's inquire what 
constraints on T uu and T vv arise when we do require S0(3) q symmetry. 



The ansatz (31) is trivially invariant under <p rotations, so the only issue is how to 
implement symmetry under the conformal Killing vector (, which in (u, v, x±, 0) coordinates 
reads 

2 , ( ® ® \ 2 2 2 \ ,9 1 + q 2 uv — q 2 x\ d 
Q = 2q x± cos (p u- — V v— + (1 + q uv + q x ± ) cos <p- sin < 



du dv J ' ' dx± x± dcj) 

(33) 

We have learned from previous sections that the useful notion of symmetry under ( involves 
(^-weights. When e has C-wdglit 4, w M has ^-weight —1, and has ^-weight —2, the 
hydrodynamical stress tensor T^ v has £- weight 2. This is not trivial to verify when the shear 
viscosity is non-zero, but it is truej^] The only sensible way to respect symmetry under £ 
prior to the collision is to demand that T^ v should again have ^-weight 2: that is, 



£ f T^ = --(V A C A )T^. (34) 



In short, the problem to be solved in this section is to find a solution to (34) of the form 



(31). For simplicity I'll also demand that the initial state preserve the Z 2 symmetry, which 



exchanges u and v: this just means that the left- and right-moving participants are identical. 

7 The Weyl covariant derivative introduced in |16j makes it substantially easier to check that has 
weight —1. 
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The answer can be anticipated from AdS/CFT: directly from [7J we can read off the result 

2g 2 £ P/ v m 2q 2 E 



T 



5(u) 



n{l + q 2 x 2 ± ) 3 
Here E is the energy in one of the shocks. 



T 

J. 71 



7i(l + q 2 xl) 3 



5(v) 



(35) 



Let me now show how (35) follows directly from (34) together with the ansatz (31). In 
(u,v,xj_,<f)) coordinates, The /x = u, v = x± component of (34) vanishes only if uT uu = 0. 
So T uu = e(x±)5(u) for some function e(x±). Plugging this expression into the /i = u, v — u 



component of (34) leads immediately to 



1 + q 2 x 2 ± )e'{x ± )S(u) + 2q 2 x ± e(x 1 _)(A5{u) + uS'(u)) = . 



(36) 



This can be simplified by using the distributional identity u5'(u) = —5{u): then (36 ) becomes 



1 + q x ± )e'(x±) + 6q x±e(x±) = 



(37) 



which can easily be solved to find the expression for T uu given in (35). A similar argument 



applies to T vv , based on the vx± and vv components of (34). It is straightforward to check 



that all other components of the equation (34) are satisfied. 



Already in [7] it was remarked that the transverse profile of T uu in ( 35 ) differs significantly 
from a highly boosted Woods-Saxon profile: See in particular figure 4 of that work. If q is 
chosen so that the energy-weighted root-mean-square (rms) value of x± is the same between 
the S'0(3)g-symmetric and Woods-Saxon profiles, then the principal difference between the 
two is that there is more weight near x± = in the SO (3) ^-symmetric profile. It is also 
significant that the large x± tail of the 50(3)g-symmetric profile is qualitatively larger than 
for Woods-Saxon: power-law falloff as compared to exponential. Ideally, one should develop 
a theory of how deviations from SO(3) q symmetry in the initial state propagate to the 
hydrodynamical stage of the collision. 



7 Plugging in approximate numbers 

Because my motivation was to understand the hydrodynamic phase of the quark-gluon 
plasma produced in heavy-ion collisions, I would now like to plug in numbers which are 
at least approximately representative of a real-world gold-gold collision at top RHIC ener- 
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gies. For some of the simpler quantities I will follow [7j: 

i = 4.3 fin £ = £ bcam = 19.7TeV /, = J- = H . (38) 
q T 4 



If the inviscid flow (22) is our goal, then it remains only to provide a value for the dimen- 



sionless quantity €q. If we want to discuss the viscous solution (26) quantitatively, then we 
must instead provide values for the dimensionless quantities T and H . 

An obvious plan for getting at e or T Q is to compute the entropy per unit rapidity in 
the fluid and then compare with phenomenological estimates of the same quantity]^] In the 
final, hadronized state, entropy is related to the number of charged tracks: 

^ re 7.5 dN *"** ^ 5000. (39) 
drj drj 



The reader interested in details of where (39 ) comes from is again referred to [7j and references 
therein, particularly [IT]. GL/V c h ar ged/c^7 is directly measurable, and for the most central 
collisions at ^/snn = 200 GeV, a reasonable figure is dN charged /dr] = 660 [IB]. This is the 
figure that went into the second approximate equality in (39)]^] 
The entropy density can be determined from the energy density: 

s = S e 3/4 where S = \f l J A = 2 A3 . (40) 

In order to compute dS/drj for the hydrodynamical flows found in previous sections, we must 
use the entropy current, 

= su» , (41) 



where s is given by (40). If M is a co-dimension 1 surface, with coordinates y a , whose 
induced metric is h a p and whose unit normal vector is everywhere timelike, then the 
entropy on M is 

S M = / d 3 y y/det hapn^. (42) 

J M 

If we take M to be a slice of constant r, and use coordinates y a = (j],x±,<f>), then the 



integral in (42) diverges because the integrand is //-independent — due precisely to the boost 
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In the interests of a compact presentation, I will not distinguish between rapidity and pseudo-rapidity. 



9 Because (39) refers to final state entropy the entropy of the fluid per unit rapidity might be significantly 
lower. Indeed, using the estimates of energy density in [TO] one finds (dS/dr])^^ ps 3000: see appendix A. 
To be more systematic in my treatment of dS/drj, I would have to introduce some definite assumptions about 
hadronization. This would take me too far afield from my main purpose, which is to obtain approximate 
numbers for the hydrodynamical flows I have found. 
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symmetry. But we can define the entropy per unit rapidity as 

/*oo /*oo 

27tt / x±dx± s T = 27rEo r / x±dx± e 3//4 w T . 
Jo Jo 



dS 
drj 



(43) 



For the inviscid flow, where we take e from (22) and u T from rt8J) and (12), the integral in 



(43) can be done explicitly, and one finds 

dS 
drj 



47r£ e 



3/4 



(44) 



There is no r dependence because entropy is conserved by inviscid flows. There is also no 



dependence on q in (44), which could be anticipated since q is dimensionful and none of the 



other quantities is. Comparing (39) and (44), and using (40), one finds 

3 ^\ 4/3 



e 



V3 \ 16 



)7T 



drj) 



880. 



(45) 



To treat viscous flow, we must decide on a value for the shear viscosity. I will take as 
a representative number the lattice result f]/s = 0.134 for SU(3) gluodynamics [20]: a bit 
larger than the value rj/s = 1/Att found in [2TJ |22] . Thus 



H = = 0.33. 

s 



(46) 



Non-zero shear viscosity implies that entropy increases with time. Because (39) refers to 



final state entropy, we should compare it with dS/ dt] on a fairly late time-slice. I will again 
evaluate entropy on a surface M at constant r. Recalling that the temperature formally 
becomes negative at large enough x±, I see that I have to cut M off at some limiting value 
x_u. To find x_u, one must solve the equation 



g* 



1 - q z T z + q z xj^ 
2qr 



(47) 



for xi», where is the unique real root of the equation T(g) = 0. In short, 

dS 



2iit / x±dx± s T 
dv Jo 



(4f 



where s T is computed as before, only using the full viscous solution (26). In order to have 
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e (GeV/fm 3 ) 



e (GeV/fm 3 ) 




Figure 1: Left: e as a function of r for different values of x± at zero viscosity, with parameters 



chosen as in (38) and (45). Right: e for non-zero viscosity, with parameters chosen as in (|38j), 
(46), and (49). The bold red line shows the dependence e = 5.4/r 3 , where e is in GeV/fm 3 



and r is in fm/c. The estimate e = 5.4 GeV/fm at r = 1 fm/c is taken from 



(48) match the result (39) with H = 0.33, one needs 



5.55. 



(49) 



This is to be compared to v^eo = 5.45 from the inviscid flow based on (45). In figure [l] 
I show the time evolution of e(r, x±) for several values of x±. For energy density I have 
used the standard units GeV/fm 3 . To convert to units of fm~ 4 , one need only recall that 
he = 0.197327 GeVfm and set h = c = 1. In figure[2]I show the direction of the hydrodynamic 
flow and contours of constant temperature for the viscous flow. In figure [3] I show the 
dependence of transverse velocity v± = tanh k and the temperature T on x± for several fixed 
values of r, using the viscous flow. 

Phenomenologically oriented readers will notice in these last plots undesirably large v± 
and T for large x±. These tails are a consequence of assuming exact conformal symmetry in 
the underlying state and of supplying initial states whose stress tensor has only power-law 
fall-off at large x±. To understand the tails better, let's consider what they look like at 



fixed r > 0. From (12) one sees immediately that v± oc l/x± for large x±. From (26) it 



follows that the temperature falls to at the limiting value x±* discussed around (48). This 



strange behavior is associated with a complete breakdown of the derivative expansion on 
which hydrodynamics is based. 

In more practical terms, the region of the flows where T < 130 MeV should probably 
be regarded as phenomenologically unuseful: 130 MeV is approximately the temperature 
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210 MeV 



130 MeV 



30 MeV 



ii (fm) 



Figure 2: The vector field w M and the temperature profile for the viscous flow, with parameters 



chosen as in (38), (46), and (49). To improve readability I have plotted not (u T ,u ) but 
instead (w r , \/ (w T ) 2 + (u- 1 ) 2 . The thick red contour is T = 130 MeV. The cooler parts of 
the plot have little to do with heavy-ion phenomenology, but they help illustrate the nature 
of the flow. 

of decoupling in a Cooper-Frye treatment, and it is significantly below T c « 170 MeV, so 
neither the conformal approximation nor the hydrodynamic approximation is any good when 
T < 130 MeV. In short, one could regard the surface T = 130 MeV as a freeze-out surface 
beyond which the true degrees of freedom are nearly free hadrons rather than a locally 
equilibrated fluid. It is worth noting that at early times, most of the energy is in the region 
where T > 130 MeV, even though the region x± < xj_* is substantially larger. In particular, 
for the viscous flow at r = 0.5fm/c, slightly more than 97% of the energy is inside the 



freeze-out surface, which is at x± = 8.4 fm, whereas T (as defined from (26)) drops to zero 
at £_u = 29.5 fm. 



The value of q in ( 38 ) is the one which makes the energy- weighted rms transverse radius 



match between the SO (3) q - invariant stress tensor profile (35) and the highly boosted Woods- 
Saxon profile of a gold nucleus. In other words, I chose the SO(3) q symmetry to be the one 
most nearly realized by the pre-collision state, and then I used that symmetry (assumed to 
be exact and to be preserved by the underlying dynamics) to constrain and after 
the collision. While this idea seems suitably straightforward for a first attempt, it could be 
in need of refinement to account for effects of evolution of the parton distribution function, 
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Figure 3: Left: The transverse velocity v± = tanh/t as a function of x± for selected values 



of r, measured in fm/c, with 1/q = 4.3 fm as in (38). Right: The local temperature T as a 
function of x± for selected values of r with non-zero viscosity, with parameters chosen as in 
(38), (46), and (49). In both plots, the red dots show where T = 130 MeV. For r ^ 4.6 fm/c, 



the temperature is below 130 MeV everywhere. 



initial state radiation, and other early time dynamics. As (44) makes clear, eo for the 



inviscid flow is entirely independent of q. Also, H is independent of q because it amounted 
simply to a choice of r\j s. Finally, Tq is only weakly dependent on q when it is determined 



by matching (48) with the observed multiplicity per unit rapidity. Thus we can think of 



varying q independently of these other parameters. An interesting quantity which depends 
strongly on q is the transverse velocity 



v± = tanh k 



2q 2 rx± 



1 + q 2 r 2 + q 2 x 



(50) 



In figure [4] I plot v± at r = 0.6 fm for several different values of q and show for comparison 
a radial flow profile studied in [2]. The comparison seems to favor 1/q somewhat larger 
than 4.3 fm. Energy density, entropy density, and temperature get smaller when 1/q is 
made larger, simply because we are expanding the flow in all space-time dimensions without 
increasing total entropy. However, the freeze-out surface changes only slowly with 1/q. It 
would be interesting to work through a more systematic study of single-particle spectra and 



Hanbury Brown-Twiss radii with (50) as an initial condition for the hydrodynamic flow: 



then one could establish a preferred value of q purely from final-state observables. 
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Xj_ (fm) 



Figure 4: The transverse velocity v± as a function of transverse radius x± at r = 0.6fm/c 
for several values of 1/g, measured in fm. The dark line is the phenomenological proposal 
v ± = tanhf^ of 0. 

8 Mapping to the covering space 

Exact solutions to the Navier-Stokes equations are few and far between. I was able to get 
hold of the one described in this paper because I imposed enough symmetry constraints so 



that there was effectively only one independent variable left, namely g defined in (22). This 
variable is essentially the only combination of the coordinates on R 3 ' 1 which is invariant under 
SO(3) q x 5*0(1, 1). It is instructive to understand what g looks like when we conformally 
embed Minkowski space in its covering space S 3 x R. An explicit mapping to accomplish 
this can be found, for example, in [23], and can be characterized as follows. Parametrize 
R 3,1 by {t,r,6,(j)), where t is the usual lab time, r = ^/(x 1 ) 2 + (x 2 ) 2 + (x 3 ) 2 , and (#,</>) 



are standard coordinates on S 2 . Parametrize S 3 x R by (cr,x,0,4>), where a is global time 
(for the R piece), (9, 0) are coordinates on S 2 as before, and \ is an angle describing the 
dimension orthogonal to the S 2 parametrized by (9, </>). Let's make all quantities on S 3 x R 
dimensionless, so that the metric is 



ds 2 = -da 2 + dx 2 + sin 2 \ (d9 2 + sin 2 9 dcj) 2 ) 



(51) 



Then the embedding of R 3,1 into S 3 x R is specified by 



qt = W sin a qr = W sin \ 



(52) 
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where 



W 



1 

cos a + cos x 



(53) 



It appears from (52) and (53) that a is an angular coordinate, to be identified modulo 2w. 
In fact, one period of a is enough to cover R 3 ' 1 . 
One can obviously embed 5 3 in R 4 as 







^sin x sin 9 cos 


<A 


2/2 




sin x sin 9 sin 





2/3 




sin x cos 9 




w 




\ cosx 


/ 



(54) 



What I have called 50 (3) g is simply the 5*0(3) which rotates yi,yz, and y± among themselves 



while leaving y 3 fixed. Slightly tedious manipulations starting from (52) and (53), together 
with x 3 = r cos 9, lead to 



2^2 



1 + q z x 



q 2 r 2 



cos a 



2qr 



^in 2 a — y\ 



(55) 



The last expression is obviously 50 (3)-invariant, since a and j/3 are separately. But 5*0(1, 1)- 



invariance, which is obvious from the middle expression in (55), is far from obvious in the 
coordinates on 5 3 x R. When r = 0, then also x± = y^ = 0, and g = cot a is a function 
simply of global time. Away from r = 0, g depends also on the "latitude" variable 2/3 on 5 3 . 
Note that slices of constant j/3 on 5 3 are two-spheres, though not the ones parametrized by 
(9,4>). The 50(3) symmetry acts to rotate these two-spheres. 

To further understand the geometry of the flow, let's start with a general description of 
conformal mappings and conformal weights. Given two manifolds with metrics, (M, g^ u ) and 
(M , g^ u ), a conformal map from M to M is a smooth bijection x^ — > x^ such that 



1 dx K dx 



A 



W 2 dx^ dx 



(56) 



where W is the conformal factor. If W — 1, then the map is an isometry. For the map (52) 



the conformal factor is given in (53). A tensor Q V a^- on M maps to a new tensor Q 
on M with conformal weight a iff 



V\_V2 — 
A t lA t 2--- 



dx Kl dx K2 dx Ul dx U2 . N N 
dx^i dxf*2 Q x Xi q x \ 2 ^«i«2- • 



(57) 
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My previous notion of weight is in a sense a special case of (|57|): Q v ^ 2 ... on R ' has £• 



weight a iff it maps to itself with conformal weight a upon the conformal map x M — » a;'* 
where d is an infinitesimal parameter which formally squares to 0. 

The four- velocity vector is naturally a tensor of conformal weight —1, because then 
its image is a timelike unit vector with respect to the new metric g^. One can check 
that for the conformal map R 3,1 ->S 3 xR discussed above, the only non-zero components 
of are in the a direction and the "latitude" direction, orthogonal to surfaces of constant 
ys ^_ One can also check that the initial shock wave states discussed in section [6]get mapped 
to shocks whose energy density is uniform across two-spheres of fixed y^ and which reach 
ys = at global time a = 0. Thus the overall picture on 5* 3 x R is as much like the original 
Bjorken picture as it can be: uniform light-like shocks collide, and a boost-invariant fluid 
results which is uniform in the transverse directions. 



9 Summary 

In order to accommodate finite transverse size and non-zero radial velocity in collisions of 
heavy ions, I propose to replace translation invariance in the transverse plane by symmetry 
under special conformal transformations, one of which is 



C - — + q 2 

dx 1 



d d 

£ijb JU _ JU JU I 



dx^ dx 1 



(58) 



where q is a parameter with dimensions of inverse length. Along with rotations around 
the beamline, these special conformal transformations fill out an 50(3) subgroup of the 
conformal group 5*0(4,2), which is an approximate symmetry of QCD at high energies. I 
denote this 5*0(3) subgroup SO(3) q . A shock- wave traveling at the speed of light in the +x 
direction which respects SO (3) q must take the form 



3 



Tun = 7Z ■ 2 2 \3 6 ( U ) ' ( 59 ) 

7r(l + q 2 x z j > d 



10 One can go further and show that if a tensor on R 3,1 respects the SO(3) q symmetry with ^-weight 
a, then mapping it with conformal weight a to S 3 x R with conformal weight a results in a tensor on 
S 3 xR which has vanishing Lie derivative under the Killing vectors that generate the SO(3) that preserves 
surfaces of constant 1/3. A special case of this general claim is that the energy density on S 3 (ignoring the 
conformal anomaly of the stress tensor) is e = W 4 e — ^-e, which is S f O(3)-invariant because both e and 
^ = 9 (sin 2 a - yl)- 1/2 are S*0(3)-invariant. 
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where E is the energy of the shock wave. Comparing with a boosted nucleus whose energy 
density is assumed to follow the Woods-Saxon profile that describes gold nuclei, one finds 
1/q — 4.3 fm in order to get the same energy-weighted root-mean-square transverse radius. 



The form (59) can be derived straightforwardly from the gauge-string duality as the dual of 



a point-sourced shock wave in AdS 5 [7] ; however, like all the results in this paper, it doesn't 
depend at all on the dynamical content of the gauge-string duality. 

The only four- velocity profile which respects SO(3) q in addition to SO(l, 1) boost invari- 
ance along the beamline and the Z 2 symmetry that reflects x 3 — > —x 3 is 



1 + q 2 r 2 + q A x 



2^2 



u 



7-1 



2gr v / T 



u 



l±v± 



qx ± 



y/l + 



(60) 



where I have written out the non-zero components of w M in the (r, rj, x±, 0) coordinate system 
and used 

1 + q 2 x\ — q 2 r 2 



9 



2qr 



(61) 



which is essentially the only combination of coordinates on R 3 ' 1 invariant under SO(3) q x 



5*0(1, 1). The four-velocity (60) might be a useful starting point for hydrodynamic simula- 
tions. In fact, a solution to the Navier-Stokes equations can be found in closed form based 



on the four- velocity (60), provided the hydrodynamic stress tensor satisfies the constraints 



of conformal invariance. Those constraints are p = e/3, vanishing bulk viscosity, and shear 
viscosity given by rj = Hoe 3//4 for some dimensionless constant Ho- In this solution, the 
temperature in the local rest frame of the fluid is 



T 



H £ 



i-(i+^Vxg,i;|;V) 



(62) 



where T is a dimensionless integration constant. Semi-realistic numbers for a central gold- 
gold collision at a/snn = 200 GeV are T = 5.55 and H = 0.33, if we choose 1/q = 4.3 fm. A 
visual presentation of the flow can be found in figure [2] the behavior of the energy density 
as a function of r for several fixed values of x± can be found in figure [TJ and plots of the 
transverse velocity v± for several different values of q can be found in figure |4j 

Although my use of the SO(3) q symmetry was inspired by head-on collisions of point- 
sourced shocks in AdS^ [7j, I do not rely on strong coupling dynamics, which has been 
argued to lead to rapidity-dependent final states [2U [25]. Instead, outside the context 
of the gauge-string duality, I am studying what might be termed "conformal collisions." 
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Suppose we have two shock waves of the form (59) colliding head-on in flat four-dimensional 
Minkowski space, and suppose the underlying dynamics of the collision is a conformal field 
theory — but let's not make any supposition about the strength of the interactions other 
than to assume that at some point after the collision, hydrodynamics applies. Finally, let's 
assume, as Bjorken suggested, that an approximate boost symmetry along the beam axis 
arises near mid-rapidity. Then, without any further knowledge of the dynamics prior to local 
equilibration, we can confidently assert that the local four-velocity in the hydrodynamical 



phase is (60) (near mid-rapidity, of course), and that the local temperature is (62). Naturally, 
my interest in these conformal collisions stems from the hope that central heavy-ion collisions 
might in some approximate sense respect SO(3) q symmetry, at least in early stages where 
the energy density is high enough for conformal invariance to be a good symmetry of the 
underlying QCD dynamics. 

The SO(3) q symmetry becomes more transparent when one maps R 3 ' 1 conformally to 
its covering space S 3 x R: it is just the rotational symmetry around one particular axis 
through the S 3 . But the SO(l, 1) symmetry is more obscure in the S 3 xR description. All 
the symmetries are manifest in an AdS$ description through the gauge-string duality. As 
remarked in 0[TT], off-center collisions of the light-like shocks described in section [6] preserves 
a U(l) subgroup out of the SO(3) q symmetry group of central collisions. (In addition, some 
discrete symmetries are preserved.) If the impact parameter is in the x 2 direction, then the 
£7(1) symmetry is the one described in (58), generalizing translations in the x 1 direction. 
This U(l) symmetry, combined with the hypothesis of boost symmetry along the beam axis, 
imposes significant constraints on the flow, but it seems in the non-central case that one can 
no longer extract the local four-velocity from symmetry considerations alone. 
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A Another entropy estimate 



It is reasonable to assume that (39 ) provides an upper bound on the entropy per unit rapidity 
of the fluid. Let's consider another estimate of entropy which is more specific to the fluid. 
In [HI] (and earlier, see e.g. [25]) one can find the estimate 



e = 5.4 



GeV 



fm d 



at r = 1 fm/c . 



(63) 



A conventional assumption is that at roughly this time, the four-velocity profile of the fluid 
is u = d T (i.e. boost invariant with no transverse flow), and e(r, x±) is proportional to a 
transverse distribution of nucleons determined by the Woods-Saxon profile: 





/ dx 3 / 


l _|_ e (\/ X3+x i- R )/ a 


e(r, x±) 


I — oo 




e(r,0) 


POD 

/ d%3 

J — oo 


/ [1 + e ilx3l ~~ R)/a ] 



(64) 



where R = 6.38 fm and a = 0.535 fm. The entropy per unit rapidity (near mid-rapidity) is 

3/4 



dS /*oo POO 

— — = 2ttt / x±dx± s(t, x_i_) = 27rrS e(r, 0) 3,/4 / x±dxj 
d'l Jo " " ~ Jo 



x ± ) 
e(r,0) 



3000, (65) 



where the final number came from plugging (63), (64), and (40) into the explicit integral in 



(65). It is reassuring that the result (65) is smaller than (39). 
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